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Abstract 


We  exercised  the  History  Variable  Reactive  Bum  (HVRB)  model  in  the  CTH  code, 
simulating  a  number  of  common  explosive  sensitivity  experiments,  including  sustained-  and 
pulsed-shock  initiation,  projectile-impact  initiation,  and  detonation  failure.  Predictions  of 
sustained-shock  initiation  are  accurate  in  accordance  with  the  calibration  of  the  model. 
Predictions  of  pulsed-shock  initiation  are  only  somewhat  more  accurate  than  those  achieved 
previously  using  Forest  Fire.  Predictions  of  projectile-impact  initiation  are  much  better.  The 
results  for  bare  explosive  show  excellent  agreement  with  experiment  for  both  flat-  and  round- 
tipped  projectiles.  With  a  2-mm  tantalum  cover,  the  agreement  is  nearly  as  good.  With  a  6-mm 
tantalum  cover,  the  agreement  deteriorates  for  flat-tipped  projectiles,  while  remaining  adequate 
in  the  case  of  round-tipped  projectiles.  The  failure  radius  and  thickness  predicted  with  CTH 
exhibit  minimam  (near  which  the  solutions  are  most  stable)  as  functions  of  zone  dimension,  but 
do  not  appear  to  converge.  These  minimum  values  agree  well  with  experiment.  Although  no 
converged  solutions  were  obtained,  the  most  stable  failure  thickness  is  approximately  equal  to 
the  most  stable  failure  radius. 
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1.  Background 


A  number  of  different  explosive  initiation  models  implemented  in  continuum  niechanics 
computer  programs  (commonly  called  hydrocodes)  have  been  offered  over  the  years.  We  have  used 
the  Forest  Fire  model  (Mador  1970, 1979;  Mader  and  Forest  1976;  Lundstrom  1988)  in  the  2DE  code 
(Kershner  and  Mader  1972)  on  many  occasions  and  developed  a  model  of  our  own  that  we 
implemented  in  a  one-dimensional  code  called  ODES  (Starkenberg  1989).  Because  2DE  is  a 
relatively  old  program,  lacking  many  modem  amenities,  including  modem  material  modeling 
capabilities,  we  were  seeking  alt^atives  for  reactive  modeling  applications. 

Several  choices  fw  the  explosive-initiation  modeler  now  exist  For  example,  anodier  early  model 
called  Ignition  and  Growth  (Lee  and  Tarver  1980)  has  been  incorporated  into  DYNA.  A  modified 
version  of  Forest  Fire  is  available  in  SMERF  (Lundstrom  1988).  More  recently,  both  Forest  Fire  and 
JTF  (a  more  sophisticated  model  named  for  its  creators  [Johnson,  Tang,  and  Forest  1985])  have  bear 
inq>lemented  in  the  three-dimensional  version  of  the  MESA/PAGOS  A  code.  In  addition  to  Forest 
Fire  and  Igmtion  and  Growth,  the  CTH  code  (Hertel  et  al.  1993)  offers  the  rather  redundantly  namaH 
History  Variable  Reactive  Bum  (HVRB)  model  (Kerley  1992).  These  models  are  all  limited  in  the 
type  of  explosive  initiation  they  apply  to.  They  are  useful  for  treating  problems  of  shock-to- 
detonation  transition  (SDT),  often  called  shock  initiaticm.  Models  for  more  complex  scenarios 
involving  deflagration-to-detonation  transition  (DDT)  have  been  attempted.  Most  of  these  are 
limited  to  specific  experimental  arrangements.  One  model  attempting  to  treat  general  geometries 
was  developed  by  Baer  and  Nunziato  (1989),  but  it  exhibits  some  difficulties  (Menikoff  1996). 

The  most  up-to-date  of  these  codes  are  MESA  and  CTH.  Because  reactive  mnHp.ling  was  only 
available  in  the  three-dimensional  v^ion  of  MESA  and  there  appears  to  be  a  genoial  preference  for 
use  of  CTH  at  the  U.S.  Army  Research  Laboratory  (ARL),  we  focused  on  the  latter.  Using  CTH, 
we  performed  computations  simulating  a  variety  of  shock-initiation  experiments  and  compared  the 
computed  results  widi  those  of  the  experiments,  as  well  as  those  obtained  firom  computations  with 
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other  codes.  The  experiments  simulated  include  sustained-  and  pulsed-shock  initiation,  projectile- 
impact  initiation,  and  detonation  failure. 

2.  Discussion  of  CTH  and  the  HVRB  Model 

CTH  is  an  ongoing  project  of  Sandia  National  Laboratories  (SNL).  It  is  intended  to  provide 
capabilities  for  modeling  dynamics  of  multidimensional  systems  with  multiple  materials,  large 
deformations,  and  strong  shock  waves.  Finite-difference  analogs  of  the  Lagrangian  equations  of 
momentum  and  energy  conservation  are  employed  with  continuous  rezoning  to  construct  Eulerian 
differencing.  Shock  and  detonation  waves  are  treated  using  the  method  of  artificial  viscosity.  CTH 
makes  use  of  analytic  (Mie-Grflneisen,  JWL,  etc.)  and  tabular  (Sesame)  equations  of  state,  as  well 
as  modem  constitutive  models  (Johnson-Cook,  Zerilli-Armstrong)  including  fracture  (void 
insertion).  In  addition,  three  reacticni  and  two  porosity  models  are  included.  These  options  provide 
an  opportunity  to  treat  complex  material  behavior,  including  melting,  vaporization,  solid-phase 
transitions,  chemical  reaction,  and  electronic  excitation  and  ionization. 

The  reaction  models  are  the  Programmed  Bum  and  CJ  Volume  Bum  models  for  detonation 
propagation,  and  the  HVRB  model  for  shock  initiation.  The  Programmed  Bum  model  forces 
detonation  at  die  characteristic  propagation  velocity  through  a  specified  portion  of  the  computational 
mesh,  while  the  CJ  Volume  Bum  model  p«mits  representation  of  self-propagating  detonation.  The 
HVRB  model  (Kerley  1992)  is  designed  to  treat  the  process  of  initiation  of  detonation  in  shock- 
loaded  high  explosives. 

In  CTH,  chemical  reaction  is  described  by  a  single-reaction  progress  variable,  X,  which  is 
assumed  to  be  a  function  of  a  history  variable,  ({). 

X  =  X(<J)), 

where 
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<|)  =  |4>(p,e,p,T,...)dT, 

0 

and  CTH  integrates  $(t)  radier  than  A,(t).  Specific  forms  for  O  and  X  are  suggested  by  the  fact  that 
the  sensitivity  of  many  explosives  to  initiation  by  pulsed  shock  waves  exhibits  a  dependrace  on  the 
square  of  the  shock  pressure  multiplied  by  the  pulse  duration.  Thus, 


P  -  Pi 

k  Pr  > 


z 


» 


and 

k  =  <|)M. 


Here,  t,,,  which  is  used  to  make  4>  nondimensional,  is  taken  to  be  equal  to  1  ps,  and  M,  pi  (the 
dueshold  pressiue  for  reaction),  p^,  and  Z  are  constants  that  may  be  calibrated  to  duplicate  die  shock- 
initiation  behavior  of  each  specific  high  explosive.  Calibration  is  accomplished  with  reference  to 
simple  sensitivity  data  characterizing  an  explosive  collected  in  the  wedge  test  Data  from  such  tests 
are  generally  presented  in  the  form  of  plots  of  distance  of  run  and  time  to  detonation  as  a  function 
of  die  initial  pressure  of  the  shock  wave.  The  formin'  of  these  is  known  as  the  Pop  plot  (Ramsay  and 
Pqpolato  1965).  Calibrations  for  several  high  explosives  have  been  obtained  and  are  disfrihntftd  with 
CTH.  In  applying  HVRB,  it  is  important  to  remember  that  teal  explosives  exhibit  modes  of 
initiation  in  addition  to  shock  initiation  as  described  by  the  model. 

As  has  become  common  practice,  CTH  treats  the  reacting  explosive  as  a  two-phase  mixture  of 
reactants  and  products.  Whm  this  is  done,  available  equations  of  state  for  the  phases  may  be  used. 
However,  since  the  indq)endent  state  variables  are  known  only  for  the  mixture,  additional 
equilibrium  conditirms  are  required  to  provide  a  closed  system  of  equations  that  can  be  solved  for 
the  dependent  state  variables.  In  some  codes,  the  phases  are  assumed  to  be  in  mechanical  and 
thermal  equilibrium  (Mader  1979).  That  is,  they  are  assumed  to  have  equal  pressures  and 
temperatures.  Generally,  the  times  required  to  achieve  thermal  equilibrium  (via  conduction)  are 
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significantly  longer  than  those  required  to  achieve  mechanical  equilibrium  (via  wave  propagation). 
Therefore,  the  condition  of  mechanical  equilibrium  can  be  satisfied  more  often  than  the  condition 
of  fiiermal  equilibrium.  Because  of  this,  die  condition  of  adiabatic  reaction  has  sometimes  been  used 
to  supplant  diermal  equilibrium  (Johnson,  Tang,  and  Forest  1985;  Starkenberg  1989).  A  significant 
disadvantage  of  any  such  procedure  is  the  requirement  for  time-consuming  iterative  computations 
in  order  to  satisfy  the  .closure  conditions.  CTH,  however,  takes  neither  of  these  approaches.  Rather, 
the  pressure,  p,  and  specific  internal  energy,  e,  of  the  mixture  are  assumed,  without  physical 
justification,  to  be  related  to  die  phase  pressures  and  specific  internal  energies  by  the  CTH  Two-State 
model: 


p(p,T,X)  =  (1  -  A)p/p,T)  +  Xpp(p,T), 
and 


e(p,T,X)  =  (1  -  X)e,(p,T)  +  Xe^(p,T), 


where  the  subscripts  r  and  p  refer  to  reactants  and  products,  respectively.  Thus,  each  phase  is 
assumed  to  have  the  density  and  temperature  of  the  mixture.  While  the  latter  assumption  represents 
thermal  equilibrium,  no  iterations  are  required  since  temperature  is  an  independent  variable  in  the 
CTH  scheme.  The  mixture  pressure  and  the  reaction  variable,  A,  have  no  physical  definition.  The 
latter  is  a  generalized  variable  that  combines  both  EOS  [equation  of  state]  and  reaction  rate 
phenomena  (Kerley  1992).  Presumably,  such  lack  of  physical  fidelity  may  be  compensated  for  in 
the  calibration  process. 

During  die  course  of  this  study,  a  new  version  of  CTH  (June  1996)  was  released.  This  version 
includes  several  new  features  including  the  Forest  Fire  and  Ignition  and  Growth  initiation  models. 
Most  importantly,  we  noticed  that  somewhat  different  results  were  obtained  using  HVRB  in  the  new 
version.  While  these  differences  were  not  large,  they  were  significant  enough  to  cause  us  to  rerun 
all  the  computations  originally  made  with  the  earlier  version.  The  results  given  in  this  report  were 
obtained  using  the  June  1996  release. 
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3.  General  Features  of  the  Simulations 


Except  in  the  case  of  the  failure-thickness  computations,  the  simulations  reported  herein  are  two- 
dimensional  axisymmetric.  Insofar  as  possible,  computational  zone  dimensions  are  kept  uniform 
and  equal  in  each  coordinate  direction  in  order  to  maintain  square  zones,  hi  all  cases,  exterior 
boundary  conditions  are  of  the  “absorbing”  type.  These  boundary  conditions  permit  mass  flow  into 
and  out  of  the  mesh  and  produce  only  weak  reflections  of  incidoit  waves. 

4.  Sustaiued-Shock  Respouse  of  PBX-9404 

The  initiation  response  of  PBX-9404  has  been  studied  in  a  wide  variety  of  expoiments.  Its 
re^nse  to  sustained-shock  waves,  detomined  in  wedge  tests  (Ramsay  and  Popolato  1965),  forms 
the  basis  of  the  Forest  Fire  and  HVRB  calibrations,  and  these  models  generally  predict  sustained- 
shock  response  accurately.  We  reconfirmed  the  accuracy  of  die  HVRB  calibration  using  two- 
dimensional  axisymmetric  simulations  and  representing  sustained-shock  initiation  by  the  impact  of 
thick  copper  flyer  plates  against  the  explosive,  hnpact  velocities  of  0.55  and  2.00  mm/ps  (giving 
pressures  near  the  upper  and  lower  extremes  of  the  Pop  plot)  were  used. 

The  thicknesses  of  the  target  explosive  was  chosen  to  be  somewhat  greater  than  the  expected 
distance  of  run  to  detonation.  The  diickness  of  the  flyo*  was  chosen  to  be  half  the  target  diickness, 
which  is  sufficient  to  prevent  encroachment  of  rarefactions  firom  the  back  of  the  flyer  into  the 
reaction  zone.  The  radii  of  the  target  and  flyer  were  chosen  equal  to  the  target  thickness  in  order  to 
prevent  encroachment  of  any  reflected  waves  near  the  axis  of  symmetry  prior  to  initiation.  Two  sets 
of  computational  zone  dimensions  were  used.  One  was  chosoi  to  give  approximately  50  zones  over 
the  expected  distance  of  run  to  detonation,  and  the  other  to  give  about  100  zones.  Thus,  0.2-  and 
0.1-mm  square  zones  were  used  for  the  longer  run  distance,  and  0.02-  and  0.01-mm  square  zones 
were  used  for  the  shorter.  The  results  obtained  were  essentially  independent  of  the  zone  dimension 
and  provided  guidance  for  the  selection  of  zoning  throughout  this  study. 
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A  sequence  of  contour  plots  for  the  0.55-mni/ps  impact  is  shown  in  Figure  1.  These  give 
pressure  contours  to  the  right  of  the  axis  and  reaction  variable  contours  to  the  left  with  the  copper 
flyer  at  the  bottom.  They  indicate  that  reaction  begins  some  distance  behind  the  impact  shock  and 
overtakes  it  to  produce  transition  to  detonation.  The  curvature  of  die  leading  wave  is  associated  with 
reflections  from  the  “absorbing”  boundary.  Plots  of  pressure  versus  distance  along  the  axis  of 
symmetry  at  various  times,  shown  in  Figure  2,  indicate  that  this  delayed  reaction  produces  a  strong 
compression  or  shock  wave  that  overtakes  the  initial  shock  (still  propagating  near  its  original 
strength).  This  behavior  is  different  from  that  observed  using  Forest  Fire,  in  which  reaction  begins 
immediately  behind  the  initial  shock  wave  and  causes  it  to  grow  as  a  square  wave  during  most  of  the 
transition  process  (Starkenberg  1993).  The  behavior  predicted  using  HVRB  is  more  consistent  with 
experimental  observations  (Wackerle  et  al.  1978;  Nutt  and  Erickson  1984).  However,  experiments 
do  not  generally  indicate  such  a  strong  overtaking  compression.  Distances  of  run  and  times  to 
detonation  were  determined  from  plots  of  pressure  versus  axial  distance  at  various  times.  Because 
of  the  interval  between  plots,  these  are  accurate  to  within  about  5%.  They  are  plotted  (with  error 
bars)  as  a  function  of  the  initial  shock  pressure  in  Figures  3  and  4.  While  the  agreement  with  the 
experimental  results  of  Ramsay  and  Popolato  (1965)  is,  of  course,  very  good,  the  run  and  time  to 
detonation  at  the  lower  pressure  are  a  little  long. 

5.  Pulsed-Shock  Response  of  PBX-9404 

In  addition  to  wedge  tests,  characterizing  the  response  of  PBX-9404  to  sustained  shock  waves, 
its  response  to  pulsed  shock  waves  has  also  been  determined  in  experiments  using  thin  Mylar  flyer 
plates  (W eingart  et  al.  1980).  The  critical  impact  velocity  for  initiation  was  determined  as  a  function 
of  flyer  thickness  for  PBX-9404  charges  measuring  25.4  mm  in  diameter  and  19.1  mm  in  length, 
laterally  enclosed  in  Plexiglas.  Forest  Fire  has  proven  inaccurate  in  predicting  these  results 
(Starkenberg  1989). 

We  used  CTTH  to  simulate  these  experiments  in  axisymmetric  computations  with  0.127-  and 
1.27-nun-thick  flyers.  In  order  to  accommodate  both  flyer  thicknesses,  we  used  a  computational 
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Figure  1.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  for  the  OJS-nun/jis 
Impact  of  a  Thick  Copper  Flyer  Against  a  Bare  PBX-9404  Charge. 
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Figure  2.  Pressure  as  a  Function  of  Axial  Distance  at  Several  Times  Within  a  PBX-94(14 
Charge  Impacted  by  a  0^5-mm/ps  Copper  Flyer. 

mesh  that  is  uniform  in  the  radial  direction  and  variable  in  the  axial  direction.  The  radial  zone 
dimension  is  0.04  mm  and  the  axial  zone  dimensions  vary  between  0.02  mm  and  0.08  mm,  with  the 
finer  zones  concentrated  near  the  impact  surface.  This  gives  approximately  six  zones  across  the 
thinner  flyer.  The  thicker  flyer  was  also  simulated  with  0.04-mm  square  zones  without  altering  the 
results. 

A  typical  sequence  of  contour  plots  leading  to  initiation  is  shown  in  Hgure  5.  The  1.27-mm- 
thick  Mylar  flyer  appears  at  the  bottom  of  the  charge.  The  pressure  contours  extend  into  the 
Plexiglas  enclosure.  Reaction  begins  when  the  wave  has  traveled  about  halfway  through  the  charge, 
and  the  transition  to  detonation  is  rapid.  A  similar  sequence  of  contour  plots  at  a  slightly  slower 
impact  velocity  is  shown  in  Figure  6.  In  this  case,  although  some  reaction  occurs,  initiation  does  not 
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Run  to  Detonation  (mm) 


Figure  3.  Comparison  of  HVRB  Predictions  ^tii  Experimental  Data  for  the  Distance  of  Run 
to  Detonation  as  a  Function  of  the  Initial  Shock  Pressure  (Pop  Plot). 
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Time  to  Detonation  (ps) 


Figure  4.  Comparison  of  HVRB  Predictions  With  Experimental  Data  for  the  Time  to 
Detonation  as  a  Function  of  the  Initial  Shock  l^essure. 
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Figure  5.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  for  the  1.5-kni/s  Impact 
of  a  1.27-nim-Thick  Mylar  Flyer  Against  a  PBX>9404  Charge. 
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Figure  6.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  for  the  1.4-kni/s  Impact 
of  a  1^7>mm*Tliick  Mylar  Flyer  Against  a  PBX-9404  Charge. 


result.  The  CTH/HVRB  critical  velocity  predictions  are  compared  with  the  experiments  and  with 
Forest  Fire  predictions  in  Figure  7.  The  experimental  results  are  fit  well  by  a  straight  line  in  a  plot 
of  the  logarithm  of  the  flyer  velocity  versus  the  logarithm  of  the  flyer  thickness.  While  HVRB 
performs  somewhat  better  than  Forest  Fire  (which  grossly  overestimates  the  sensitivity  of  the 
explosive),  it  still  underestimates  the  sensitivity  considerably. 


Figure  7.  Comparison  of  HVRB  and  Forest  Fire  Predictions  With  Experimental  Data  for  the 
Initiation  of  PBX-9404  Charges  by  the  Impact  of  Thin  Mylar  Flyers. 


6.  Projectile-Impact  Initiation  of  PBX-9404 


The  initiation  of  explosives  due  to  projectile  impact  has  been  the  subject  of  extensive 
experimental  research,  hi  this  case,  initiation  is  influenced  by  rarefactions  emanating  fixim  tiie  outer 
surface  of  the  projectile  that  reduce  reaction  rates  and  attenuate  impact  shock  waves.  A  study  by 


13 


Bahl,  Vantine,  and  Weingart  (1981)  produced  critical  velocities  for  impact  of  flat-  and  round-tipped 

cylindrical  steel  projectiles  against  bare  and  tantalum-covered  PBX-9404  charges.  The  cover  plates 

♦ 

in  these  experiments  were  either  2-  or  6-mm  thick.  The  results  of  this  study  and  our  modified 
Jacobs-Roslund  fits  to  the  data  are  summarized  in  the  appendix. 

We  made  axisymmetric  simulations  of  some  of  these  experiments,  representing  the  flat-tipped 
projectiles  as  right  circular  cylinders  with  unit  length-to-diameter  ratios  and  adding  hemispheres  to 
the  leading  edges  to  produce  round-tipped  projectiles.  We  maintained  a  fixed  relationship  between 
the  projectile  and  the  target,  keeping  the  target  diameter  and  depth  equal  to  five  times  the  projectile 
diameter.  TTie  computational  zoning  was  varied  with  the  projectile  diameter.  0.1 -mm,  0.2-mm,  or 
0.5-mm  square  zones  were  selected  as  appropriate. 

In  Older  to  limit  the  number  of  computations,  we  considered  only  two  projectile  diameters  (near 
the  upper  and  lower  end  of  the  experimental  range)  for  each  of  the  three  target  configurations.  We 
varied  the  projectile-impact  velocity  in  order  to  determine  critical  values  for  initiation.  Typictil 
sequences  of  contour  plots  are  shown  in  Figures  8-15. 

Figure  8  shows  a  sequence  of  pressure  and  reaction  variable  contour  plots  for  initiation  of  a  bare 
PBX-9404  target  by  the  impact  of  a  1 6-mm-diameter  flat-tipped  projectile  at  0.60  mm/p  s.  This  near- 
critical  impact  produces  an  iiutial  shock  pressure  of  about  3.9  GPa.  hi  this  case,  the  impact  shock 
propagates  approximately  8.3  mm  into  the  target  as  a  planar  wave  (i.e.,  without  attenuation  at  the 
axis)  before  transitioning  to  detonation.  This  distance  of  run  to  detonation  is  less  than  10%  longer 
than  that  associated  with  the  computed  Pop  plot,  indicating  minimal  influence  of  the  rarefaction  on 
the  reaction  zone.  When  the  impact  velocity  of  the  same  projectile  is  reduced  to  0.55  mm/ps, 
detonation  fails  to  initiate  in  the  same  target,  as  shown  in  the  sequence  of  pressure  and  reaction 
variable  contour  plots  of  Figure  9. 

Results  obtained  with  the  projectile  diameter  reduced  to  2  mm  are  shown  in  Figure  10.  In  this 
case,  a  near-critical  1.5-mm/ps  impact  produces  an  initial  shock  pressure  of  about  13.7  GPA,  leading 
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Figure  8.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  Initiation  of  a 
Bare  PBX-9404  Target  by  the  Impact  of  a  16*nun-Diameter  Flat*Tipped  Steel 
Projectile  at  0.60  mm/ps. 
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Figure  9.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  Initiation 
Failure  in  a  Bare  PBX>9404  Target  Following  Impact  by  a  16-mm-Diameter  Flat- 
Tipped  Steel  Projectile  at  0^5  nun/ps. 
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Figure  10.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  Initiation  of 
a  Bare  PBX-9404  Target  by  tiie  Impact  of  a  2-nun<Diameter  Flat>Tipped  Steel 
.  Projectile  at  1.5  mm/us. 
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Figure  11.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  Initiation  of 
a  Bare  PBX>9404  Target  by  the  Impact  of  an  18>mm-Diameter  Round-Tipped  Steel 
Projectile  at  1.0  nun/ps. 
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Figure  12.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  the  Initiation 
of  a  PBX-9404  Target  With  a  2-nini-Thick  Tantalum  Cover  by  the  Impact  of  a 
20-mm<Diameter  F1at*Tipped  Steel  Projectile  at  0.7  mm/ps. 
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Figure  13.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  the  Initiation 
of  a  PBX-9404  Target  With  a  2-nun-Thick  Tantalum  Cover  by  the  Impact  of  a 
4-mm-Diameter  Flat-Tipped  Steel  Projectile  at  2.0  mm/jis. 
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Figure  14.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  the  Initiation 
of  a  PBX*9404  Target  With  a  6-nun-Tliick  Tantalum  Cover  by  the  Impact  of  an 
18>nim>Diameter  Flat-Tipped  Steel  Projectile  at  1.0  mm/ps. 
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Figure  15.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  Showing  the  Initiation 
of  a  PBX-94()4  Target  With  a  6>nun-Thick  Tantalum  Cover  by  the  Impact  of  a 
16-mm«Diameter  Round-Tipped  Steel  Projectile  at  2.0  mm/^s. 
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to  detonation.  The  run  to  detonation  is  about  1.5  mm,  or  nearly  50%  greater  than  that  associated 
with  the  Pop  plot  The  influence  of  the  rarefaction  on  the  reaction  zone  is  significant 

With  round-tipped  projectiles,  rarefactions  enter  the  reaction  zone  immediately  after  impact 
leading  to  hi^er  critical  velocities.  Figure  1 1  shows  a  sequence  of  pressure  and  reaction  variable 
contour  plots  for  initiation  of  a  bare  target  by  the  impact  of  an  18-mm-diameter  flat-tipped  projectile 
at  1.0  mm/ps.  The  initial  pressure  is  about  7.8  GPa,  and  the  run  to  detonation  of  about  6.0  mm  is 
more  than  twice  the  Pop  plot  value. 

The  effect  of  a  cover  plate  is  to  allow  attenuation  of  the  impact  shock  wave  by  the  overtaking 
rarefaction  prior  to  its  entry  into  the  explosive  target,  again,  leading  to  higher  critical  velocities.  If 
the  ratio  of  the  cover  thickness  to  the  projectile  diameter  is  small  enough,  the  central  portion  of  the 
shock  wave  produced  by  the  impact  of  a  flat-tipped  projectile  will  remain  planar  and  unattenuated 
as  it  enters  the  explosive.  This  is  the  case  for  a  20-mm-diameter  flat-tipped  projectile  and  a 
2-mm-fhick  tantalum  cover.  Results  obtained  witii  a  projectile  velocity  of  0.7  mm/ps  are  shown  in 
Figure  12.  hiitiation  occurs  within  the  planar  region  with  a  run  to  detonation  comparable  to  that  of 
the  Pop  plot. 

With  a  small  diameter  projectile,  or  thick  cover  plate,  the  planar  region  of  the  impact  ^ock  may 
be  totally  attenuated  within  die  cover.  For  example,  when  the  projectile  diameter  is  reduced  to  4  mm 
with  the  same  cover  plate,  an  impact  velocity  of  2.0  mm/ps  or  greater  is  required  to  produce 
initiation.  Results  are  shown  in  Figure  13.  In  this  case,  &e  run  to  detonation  is  longer  than  that 
associated  with  the  Pop  plot  Similar  results  are  obtained  in  the  case  of  a  18-mm-diameter  flat- 
tipped  projectile  and  a  6-mm-thick  tantalum  cover.  With  an  impact  velocity  of  1.0  mm/ps,  most  of 
the  planar  region  of  the  impact  shock  is  attenuated  within  the  cover  as  shown  in  Figure  14. 

When  the  effects  of  the  cover  plate  are  combined  with  those  of  a  round  tip,  even  higher  projectile 
velocities  are  required  to  produce  detonation.  Results  for  a  16-mm-diameter  projectile  with  a 
velocity  of  2.0  mm/ps  impacting  a  6-mm-thick  cover  plate  are  shown  in  Figure  15. 


23 


Critical  velocities  are  plotted  as  functions  of  projectile  diameter  in  Figures  16-20.  The  fits  to 
experimental  data  are  shown  as  solid  lines  bounded  by  fiUed  symbols,  while  the  computational 
results  arc  plotted  as  opoi  symbols  connected  by  dashed  lines.  The  results  for  bare  explosive  show 
excellent  agreement  with  experiment  for  both  flat-  and  round-tipped  projectiles,  as  illustrated  in 
Figure  16.  With  a  2-nim  tantalum  cover,  the  agreement  is  nearly  as  good,  as  shown  in  Figure  17. 
Wifli  a  6-mm  tantalum  cover,  the  agreement  deteriorates  for  flat-tipped  projectiles  while  remaining 
adequate  in  the  case  of  round-tipped  projectiles,  as  shown  in  Figure  18.  The  same  data  are  plotted 
again  in  Figures  19  and  20,  which  show  the  trends  with  increasing  target  cover  for  each  projectile 
tip  shape. 


Figure  16.  Comparison  of  HVRB  Predictions  With  Experimental  Data  for  the  Critical  Impact 
Velocity  as  a  Function  of  the  Diameter  of  Cylindrical  Steel  Projectiles  Against  Bare 
PBX-9404  Charges. 
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Figure  17.  Comparison  of  HVRB  Predictions  With  Experimental  Data  for  the  Critical  Impact 
Velocity  as  a  Function  of  the  Diameter  of  Cylindrical  Steel  Projectiles  Against 
2>nmi-Thick  Tantalum-Covered  PBX-9404  Charges. 


Figure  18.  Comparison  of  HVRB  Predictions  With  Experimental  Data  for  the  Critical  Impact 
Velocity  as  a  Function  of  the  Diameter  of  Cylindrical  Steel  Projectiles  Against 
6-mm-Thick  Tantalum-Covered  PBX-9404  Charges. 
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Projectile  Velocity  (mm/MS)  ij  Projectile  Veiccity  (mnVus) 


Comparison  of  HVRB  Predictions  With  Experimental  Data  for  the  Critical  Tmpapt 
Velodty  as  a  Function  of  the  Diameter  of  Flat-Tipped  Cylindrical  Steel  Projectile 
Against  Bare  and  Tantalum-Covered  PBX-9404  Charges. 


Figure  20.  Comparison  of  HVSB  Predictions  With  Experimental  Data  for  the  Critical  Impact 
Velocity  as  a  Function  of  the  Diameter  of  Round-Tipped  Cylindrical  Sted 
Projecffles  Against  Bare  and  Tantalum-Covered  PBX-9404  Charges. 
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7.  Detonation  Failure  in  Composition  B  Charges 


It  is  well  known  that  explosive  charges  can  only  sustain  detonation  when  their  lateral  dimensions 
are  sufficiently  large.  When  this  condition  is  not  met,  incipient  detonations  fail  as  a  result  of  the 
effects  of  rarefactions  that  encroach  upon  their  reaction  zones.  Thus,  for  a  particular  explosive,  a 
failure  diameter  (or  radius)  may  be  determined  for  cylindrical  charges,  and  a  failure  thickness  may 
be  determined  for  laminar  charges.  The  phenomenon,  which  is  of  interest  in  its  own  right,  has 
consequences  for  initiation  of  explosives  by  small  projectiles  such  as  shaped-charge  jets  (Chick  et  al. 
1985,  1989;  Starkenberg,  Lawrence,  and  Dorsey  1994)  and  affects  the  ability  of  propagating 
detonations  to  turn  comers  in  explosive  charges. 

Experimental  data  describing  the  failure  diameter  (Campbell  and  Engelke  1976;  Dobratz  1985) 
and  thickness  (Urizar  1978;  Gibbs  and  Popolato  1980)  of  several  explosives  ate  presented  in  Table  1. 
The  last  column  of  the  table  gives  the  ratio  of  the  failure  radius  to  the  failure  thickness.  While  (he 
ratio  varies  between  0.64  and  1 .21,  in  most  cases,  it  lies  very  close  to  unity.  Since  these  experimrats 
were  conducted  by  different  researchers  at  different  times,  the  explosives  may  not  be  identical,  hi 
fact,  only  Composition  A-3  (for  which  a  ratio  of  0.96  applies)  appears  to  have  the  same  formulation 
and  density  in  both  experiments.  Thus,  the  data  are  generally  consistent  witih  the  hypothesis  the 
failure  thickness  and  radius  are  equal. 


Table  1.  Experimental  Failure  Thickness  and  Radius  for  Selected  Explosives 


Explosive 

Density 

(g/cm^) 

Failure 

Thickness 

(mm) 

Explosive 

Density 

(g/cm^) 

Failure 

Radius 

(mm) 

Radius/ 

Thickness 

(mm/mm) 

PBX9404 

1.83 

0.92 

PBX9404 

0.59 

Comp.  A-3 

1.63 

1.14 

Comp.  A-3 

1.10 

Cyclotol  75/25 

1.75 

3.02 

Cyclotol  77/23 

1.74 

3.00 

Comp.  B-3 

1.72 

1.88 

Comp.  B-3 

— 

2.00 

1.06 

Octol  75/25 

1.79 

2.86 

Octol  75/25 

1.81 

3.20 

1.12 

Pentolite 

1.70 

2.78 

Pentolite 

— 

3.35 

1.21 
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The  reaction  zone  for  fully  developed  detonations  is  on  the  order  of  a  millimeter  for 
Composition  B  and  is  not  well  resolved  in  most  computations.  Such  lack  of  computational 
resolution  affects  the  accuracy  of  detonation  failure  predictions.  However,  sufficiently  fine  zoning 
should  produce  accurate  solutions  that  are  independent  of  zone  size.  Mader  (1979)  has  reported 
accurate  predictions  of  failure  diameter  in  several  explosives  using  the  Forest  Fire  model  in  the  2DL 
and  2DE  codes.  Lundstrom  (1993)  has  produced  less  successful  predictions  of  the  failure  diameter 
of  Composition  B  using  his  modified  version  of  Forest  Fire  in  the  SMERF  code.  He  observed  rapid 
increases  in  failure  diameter  with  zone  size  at  the  smallest  zone  sizes  he  used.  More  recently,  we 
obtained  accurate  predictions  of  the  failure  radius  and  thickness  of  Composition  B  that  are 
independent  of  computational  zone  size  using  2DE  and  the  original  version  of  Forest  Fire 
(Starkenberg  and  Dorsey  1994).  These  indicate  that  the  failure  thickness  of  an  explosive  is  generally 
equal  to  its  failure  radius. 

In  order  to  obtain  an  understanding  of  the  utility  of  HVRB  in  predicting  detonation  failure,  we 
made  both  axisymmetric  and  plane-strain  computations  to  predict  the  failure  radius  and  thickness 
of  Composition  B.  The  simulation  consists  of  Composition  B  booster  and  main  charges  surrounded 
by  air.  Initiation  is  achieved  by  means  of  programmed  bum  detonation  in  the  booster  charge. 
Reaction  in  the  main  charge  is  described  using  HVRB.  We  made  use  of  the  CTH  velocity  addition 
feature  to  subtract  an  axial  velocity  approximately  equal  to  the  detonation  velocity  throughout  the 
mesh  at  a  specified  time.  This  serves  to  freeze  a  detonation  after  it  has  propagated  a  desired  distance 
(in  our  case,  approximately  halfway  through  the  mesh).  When  the  detonation  fails,  the  decaying 
wave  is  swept  off  the  bottom  of  the  mesh.  We  modified  an  input  template  suggested  by  Kerley 
(1995)  for  the  purpose  of  initializing  the  computations.  As  the  charge  radius  or  thickness  is  varied, 
the  dimensions  of  the  computational  mesh,  the  length  of  the  booster  charge,  and  the  time  for  the 
velocity  addition  are  varied  proportionally  with  a  fixed  zone  size. 

Typical  sequences  of  contour  plots  for  failure-radius  computations  are  shown  in  Figures  21  and 
22.  For  these  computations,  we  used  zone  sizes  ranging  from  0.04-  to  0.20-mm  square.  Figure  21 
shows  detonation  propagation  through  7.0  ps  in  a  2.04-mm  radius  charge  using  a  zone  size  of 
0.12  mm.  A  velocity  of  7.65  mm/ps  is  subtracted  from  the  grid  at  1.1  ps.  After  this,  the  detonation 
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Figure  21.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  for  Detonation 
.Propagation  in  a  2.04-nun-Radius  Cylindrical  Composition  B  Charge. 
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Figure  22.  Sequence  of  Pressure  and  Reaction  Variable  Contour  Plots  for  Detonation  Failure 
in  a  1.92>mm-Radius  Cylindrical  Composition  B  Charge. 


front  drifts  slowly  downward  in  the  grid  because  its  propagation  velocity  is  a  little  lower  than  the 
subtracted  velocity.  In  a  similar  computation  with  a  1.92-mm  radius  charge,  the  detonation  failed 
after  about  1.5  ps,  as  shown  in  Figure  22.  All  activity  is  swept  off  the  mesh  by  3.0  ps. 

We  have  plotted  the  CTH  failure-radius  predictions,  along  with  the  results  we  obtained  with 
2DE,  as  well  as  those  of  Lundstrom,  as  functions  of  zone  dimension  in  Figure  23.  While  our  2DE 
computations  converged  nicely  with  decreasing  zone  size,  Lundstrom  did  not  observe  such 
convergence  in  his  computations,  which  exhibit  a  minimum  in  predicted  failure  radius.  CTH  also 
produces  solutions  that  appear  to  diverge  with  decreasing  zone  size,  although  not  as  rapidly  as 
Lundstrom’s.  The  apparent  divergence  does  not  imply  that  the  solutions  never  converge,  only  that 
they  do  not  converge  for  zone  sizes  that  are  practical  to  use.  The  failure  radius  predicted  with  CTH 
also  exhibits  a  minimum,  near  which  the  solution  is  most  stable  (least  sensitive  to  variations  in  zone 
size),  at  a  zone  dimension  of  0.12  mm.  This  minimum  value  agrees  well  with  the  experimental 
failure  radius  (Campbell  and  Engelke  1976;  Dobratz  1985).  The  CTH  most  stable  zone  size  is 
somewhat  larger  than  that  determined  by  Lundstrom  and  the  zone  size  required  for  convergence  with 
2DE.  The  most  stable  zone  size  can  be  applied  to  other  configurations  in  which  detonation  failure 
is  an  issue  (in  particular,  to  shaped-charge  jet  attack). 

For  die  failing  thickness  computations,  we  used  zone  sizes  ranging  from  0.01-  to  0.20-mm 
square.  Results  obtained  are  plotted  in  Figure  24.  In  this  case,  the  most  stable  zone  size  is  about 
0.02  mm,  much  smaller  than  in  the  failure-radius  case.  The  associated  failure-thickness  prediction 
is  a  little  larger  than  the  experimental  value  (Urizar  1978;  Gibbs  and  Popolato  1980).  However,  the 
experimental  method  is  indirect  and  renders  the  latter  value  suspect.  Agreement  with  the 
experimental  failure  radius  (which  may  more  accurately  represent  the  failure  thickness)  is  better. 

8.  Summary  and  Conclusions 

We  exercised  die  HVRB  model  in  the  CTH  code,  simulating  a  number  of  common  explosive 
sensitivity  experiments,  including  sustained-  and  pulsed-shock  initiation,  projectile-impact  initiation, 
and  detonation  failure. 
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Failure  Thickness  (mm) 


Zone  Dimension  (mm) 


Figure  23.  Comparison  of  HVRB  and  Forest  Fire  Predictions  for  Composition  B  Failure 
Radius  as  a  Function  of  the  Computational  Zone  Size. 


Figure  24.  HVRB  Predictions  for  Composition  B  Failure  Thickness  as  a  Function  of  the 
Computational  Zone  Size. 
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Predictions  of  sustained-shock  initiation  are  accurate  in  accordance  with  the  calibration  of  the 
model.  Predictions  of  pulsed-shock  initiation  are  only  somewhat  more  accurate  than  those  achieved 
previously  using  Forest  Fire. 

Predictions  of  projectile-impact  initiation  are  much  better.  The  results  for  bare  explosive  show 
excellent  agreement  with  experiment  for  both  flat-  and  round-tipped  projectiles.  With  a  2-mm 
tantalum  cover,  the  agreement  is  nearly  as  good.  With  a  6-mm  tantalum  cover,  the  agreement 
deteriorates  for  flat-tipped  projectiles,  while  remaining  adequate  in  the  case  of  round-tipped 
projectiles. 

The  failure  radius  predicted  with  CTH  exhibits  a  minimum  (near  which  the  solution  is  most 
stable)  as  a  function  of  zone  dimension,  but  does  not  appear  to  converge.  This  minimum  value 
agrees  well  with  experiment  Results  for  failure  thickness  are  similar.  The  most  stable  value  is  a 
little  larger  than  the  experimental  failure  thickness  and  roughly  equal  to  the  experimental  failure 
radius. 

Difficulties  associated  with  resolution  of  the  reaction  zone  may  limit  the  accuracy  of  predictions 
in  problems  where  large  amounts  of  reaction  occur  over  short  distances.  Cognizance  of  these 
difficulties  is  important  in  conducting  simulations  of  detonation  fulure,  comer  turning,  and  initiation 
by  projectiles  having  diameters  smalls  than  the  explosive’s  failure  diameter. 
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Appendix: 

Modified  Jacobs-Roslund  Fits 
to  PBX-9404  Projectile-Impact  Initiation  Data 
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The  widely  used  Jacobs-Roslund  fonnula  giving  the  critical  impact  velocity  for  initiation  of  an 
explosive  charge  as  a  function  of  projectile  diameter  may  be  written 


Vcrit  = 


AdlB) 

d'/2  ^  d, 


where  d  is  the  projectile  diameter,  h  is  die  thickness  of  the  plate  covering  the  explosive  charge;  and 
A,  B,  and  C  are  constants  used  to  rit  the  equation  to  available  data.  In  fact,  there  are  only  two 
indqiendent  constants,  A(1  +  B)  and  C.  The  constant  B,  which  vanishes  for  flat-tipped  projectiles, 
is  a  convenience  used  to  account  for  the  effect  of  projectile  tip  shape,  llie  constant  A  is  dien 
calibrated  using  data  for  flat-tipped  projectiles  and  bare  explosive  charges.  Hie  constant  C  is  clearly 
associated  with  the  effect  of  die  cover  plate. 

B^,  Vantine,  and  Weingart‘  produced  an  excellent  set  of  data  idmtif^g  initiation  thresholds 
for  bare  and  tantalum-covered  PBX-9404  charges  attacked  by  flat-  and  round-tipped  projectiles. 
Liddiard  and  Roslund*  have  summarized  the  data  (including  that  for  PBX-9404)  and  provided  fit 
constants  for  a  number  of  explosives.  Starting  with  the  bare-charge  data  for  flat-tipped  projectiles, 
they  identified  pairs  of  data  points  near  tiie  tiireshold  that  include  one  go  and  one  no-go.  They  then 
det^mined  the  midpoint  of  each  pair  by  averaging  the  associated  impact  velocities  and  the  inverse 
square  roots  of  the  associated  projectile  diameters.  They  then  computed  the  constant  A  for  each  pair. 
Since  the  standard  deviation  was  small,  they  used  the  average  value  of  A  for  the  fit  Following  a 
similar  procedure  using  the  round-tipped  projectile  data,  they  detomined  values  for  the  constant  B. 
Since  B  appeared  to  vary  witii  projectile  diamet^:,  at  least  at  larger  diameters,  they  argued  that  this 
reflected  a  ch^ge  size  effect  in  the  experiments  and  that  B  should  be  determined  fi-om  data  for 
projectiles  having  diameters  of  7  mm  or  less.  Similar  determinations  of  C,  however,  indicated  tiiat 


^  Bahl,  K.  L.,  H.  C.  Vantine,  and  R.  C.  Weingart  “The  Shock  Initiation  of  Bate  and  Coveted  Explosives  by  Projectile 
Ittq)act”  Proceedings  of  the  Ninth  Symposium  (International)  on  Detonation,  pp.  604-617,  August  1985. 

^  Iiddiatd,T.P.,andL.A.Roslund.  “Ptojectile/FragmentInq)act  Sensitivity  of  Explosives.”  NSWC-IR89-184,  Naval 
Sutface  Weapons  Center,  1993. 
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it  depends  on  tip  shape  and  is  not  constant.  The  midpoint  data  used  in  this  procedure  are  given  in 
Tables  A-1  and  A-2. 


Table  A-1.  Critical  Velocities  for  Impact  of  Flat-Tipped  Projectiles 


1  Bare  Charge 

2-mm  Tantalum  Cover 

6-mm  Tantalum  Cover  || 

d 

^crit 

d 

Vcrit 

d 

Vcri.  1 

(mm) 

(mm/ps) 

(mm) 

(mm/ps) 

(mm) 

(mm/ps)  1 

1.13 

2.070 

3.684 

MEM 

1.28 

1.910 

3.724 

■H 

MEM 

1.49 

1.620 

3.913 

1.740 

10.00 

1.26  1 

2.20 

1.380 

3.968 

1.840 

15.01 

0.95 

3.14 

1.130 

4.675 

1.605 

18.00 

0.84 

3.68 

1.085 

5.656 

1.420 

— 

— 

7.73 

0.745 

7.631 

1.130 

— 

— 

8.70 

0.695 

8.858 

1.040 

— 

— 

1  12.76 

0.575 

12.986 

0.840 

— 

— 

1  17.73 

0.500 

20.291 

0.650 

— 

— 

Table  A-2.  Critical  Velocities  for  Impact  of  Round-Tipped  Projectiles 


1  Bare  Charge 

2-mm  Tantalum  Cover 

— _ _ _ 

6-mm  Tantalum  Cover  | 

1  ^ 

Vcrit 

d 

VcH, 

d 

Vcrit 

(mm) 

(mm/ps) 

(mm) 

(mm/ps) 

(mm) 

(mm/ps) 

3.968 

1.910 

10.0 

15.0 

Msm< 

4.234 

1.785 

15.0 

Bl 

20.0 

■Hi 

4.451 

1.845 

17.3 

mem 

23.5 

1.50 

4.747 

1.720 

20.0 

1.40 

25.0 

1.44 

5.373 

1.615 

26.0 

1.25 

31.0 

1.23 

6.575 

1.530 

— 

— 

— 

— 

10.750 

1.210 

— 

— 

— 

— 

11.413 

1.185 

— 

— 

— 

— 

12.015 

1.160 

— 

— 

— 

— 

14.080 

1.100 

— 

— 

— 

_ 

16.458 

1.040 

— 

— 

— 

_  1 

18.904 

0.960 

— 

— 

— 

—  1 
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Using  all  of  this  data,  we  found  that  better  fits  could  be  obtained  in  the  form 


If  a  =  0.5  and  b  =  1.0,  this  is  identical  to  the  Jacobs-Roslund  form.  The  constants  we  used,  along 
with  corresponding  Jacobs-Roslund  constants,  are  shown  in  Table  A-3. 


Table  A-3.  Critical  Velocity  Fitting  Parameters 


A 

(mm^/us) 

a 

B 

b 

Present 

2.109 

0.51 

1.098 

0.45 

Flat 

J-R 

2.09 

0.50 

1.90 

1.00 

Present 

3.338 

0.42 

1.156 

0.40 

Round 

J-R 

3.78 

0.50 

0.630 

1.00 
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